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Summary 

During this giant period we have investigated two different aspects of active 
antenna array design. The first of these involved the study of transition between 
monolithic microwave integrated circuits and rectangular waveguides. The second 
aspect of the study was concerned with the investigation of crosstalk in 
multiconductor transmission lines. 

We have formulated the boundary value problem associated with a 
discontinuity in a microstrip line. This entailed, as a first step, the derivation of the 
propagating as well as evanescent modes of a microstrip line. We have completed this 
task and have also derived the solution to a simple discontinuity problem: change in 
width of the center strip. Currently, we are in the process of refining this solution. 

As for the multiconductor transmission line problem, we have developed a 
computer algorithm for computing the crosstalk noise from the signal to the sense lines. 
The computation is based on the assumption that these lines are terminated in passive 
loads. During the next grant period we plan to generalize this program to the case 
where the lines are terminated by active logic gates. 
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1. CHARACTERIZATION OF MMIC DEVICES FOR ACTIVE ARRAY ANTENNAS 
1.1. Introduction 

During the past year, the University of Illin ois has been involved in the study 
of transitions between monolithic microwave integrated circuits (MMICs) and 
rectangular waveguide. This problem occurs in the design of active array antennas, 
where the gain and amplitude of the individual array elements are controlled by 
MMIC’s. In this report, we will su mm arize the progress that has been made on this 
problem to-date. 

Transitions between MMICs and rectangular waveguides should have several 
properties. First, they should have a low insertion loss. Next, they should be 
structurally and thermally sound. Finally, these transitions should maintain the 
inherent reproducibility of the integrated circuit. A number of possible candidates for 
these transitions are shown in Figures 1 and 2. In order to determine which of these 
various transitions are likely to offer the best characteristics, it is necessary first to 
analyze them for insertion loss and return loss. We propose to achieve this by first 
breaking up a transition into a succession of abrupt discontinuities. If we can calculate 
the S-parameters of a single discontinuity, we should then be able to determine the 
characteristics of a cascade of such discontinuities. The purpose of this study will 
initially be to develop the tools necessary to study these abrupt discontinuities in 
printed circuits. 

Abrupt printed circuit discontinuities may take a number of forms. The 
theory which is developed is general for many different types of discontinuities, as will 
be discussed later. For the present, however, it is necessary to begin by studying the 
most simple example of this discontinuity one can find. This turns out to be an abrupt 
change in the strip width of a shielded microstrip. A cross section of uniform shielded 
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Figure 1. Typical substrate-to- waveguide transitions. 
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Figure 2. Transition between fin-line and rectangular waveguide. 
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microstrip is shown in Figure 3, and a discontinuity in the strip width is shown in 
Figure 4. By studying this structure, it is hoped that the tools necessary to study a 
more realistic transition, such as the fin-line stepped taper shown previously in 
Figure 2, will be developed. 

The method to be used in the analysis of abrupt discontinuities involves mode 
matching in the plane of the discontinuity. In order to achieve this, one begins by 
generating the dominant and first few higher-order modes in each of the two microstrip 
lines. Next, one matches the tangential electric and magnetic field components in the 
plane of the discontinuity. Finally, one calculates the mode coefficients, which yield S- 
parameters and equivalent circuits of the discontinuity. 

The study begins now r with the analysis of the uniform microstrip. 


1.2. Analysis of a Uniform Microstrip 

The first step in the analysis of a uniform microstrip involves finding a Green’s 
function which relates currents on the strip to the electric fields at all other points in 
the cross section of the line. This is accomplished by solving a two dimensional 
Helmholtz equation in regions I and II 


(V 2 + *, 2 ) 


= 0 


*. = e, 


€: = 


e, x >0 
e,e r x <0 


(la) 

(l.b) 


(lx) 


In the above equations, i denotes the region, and 0; and are the electric and magnetic 
scalar potentials in each region. 

It is convenient to solve the above equation in the spectral domain, so one takes 
the Fourier transform of all potentials and fields as 



6 


I 

h 

i 


i 

d 

T 



PEC SHIELD 
REGION I 


2b 


REGION IE 


Figure 3. Cross section of shielded microstrip. 



Figure 4. Top view of a discontinuity in the strip width on a microstrip. 
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By matching boundary conditions on the shield walls and on the center conductor, one 
obtains an equation of the form 


Z 22 (<*„,/ 3) Z„(a„,J 3) 




E z (a n ,/3) 

Z«(a n ,j3) Z xx (a„ ,0) 


Jx (a n ) 


E x (a n ,P) 


This equation relates the current in the plane y = 0 to the electric fields in this plane. 
The dyadic Green’s function, Z^( <*„ ,/3), consists of relatively simple expressions in 
terms of hyperbolic sines and cosines. 

Before proceeding to the solution of this equation for 0, one should first look at 
the methods for generalizing the above Green’s function. The above method is 
satisfactory for configurations with one dielectric layer and one strip. But a method 
developed by T. Itoh [l] gives a technique for generating a Green’s function for 
structures with an arbitrary number of dielectric layers and conducting strips. His 
technique is called the Spectral Immitance Approach. While a detailed explanation of 
the technique is beyond the scope of this paper, it may be stated that the the technique 
involves the separation of the fields in each dielectric region into TE y and TM y 
components in the spectral domain, and the subsequent formulation of analogous 
transmission lines in the y direction for the decoupled TE and TM fields. This 
generalization extends the applicability of the techniques presented in this paper to a 
large class of printed circuits. 

Next, one must find the solution to the integral equation shown in Equation (3). 
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The solution is effected by the spectral Galerkin technique [2\ in which a moment 
method solution is brought about in the spectral domain. Hence, the current on the 
strip is expanded in terms of basis functions, which are non-zero only on the strip, as 


JXx ) 

J x (x) 


E 

i=l 


E 

i =1 


COS ( (t— 1) 77 (xA) — l) 

Vl —(x/s) 2 

, sin (i it (x /s ) — l) 

' Vl - (x/s ) 2 


(4a) 

(4.b) 


One then takes the Fourier transform of these currents and substitutes the result into 
Equation (3). Next, the inner product of the resulting equation is taken with the 
individual basis functions, and Parseval’s theorem is used to eliminate the right-hand 
side. The resulting matrix equation may be solved by setting the determinant equal to 
zero and solving for j8 with Newton’s method. 

The accuracy of this method is determined by the number of basis functions 
and the number of spectral terms one can calculate within a reasonable amount of 
computer time. For the dominant mode and lower-order evanescent modes, a small 
number of basis functions and spectral terms is probably satisfactory. For evanescent 
modes of larger order, it is likely that more basis functions and spectral terms are 
required to achieve reasonable accuracy. 

Sample calculations were carried out in order to calculate the dominant and 
first two evanescent modes of microstrip. These results are shown in Figure 5. One 
would like to compare these modes to other results, but little data are available on 
microstrip modes. The microstrip mode calculations that are available deal with 
propagating modes, but do not present data on evanescent modes [3l[4]. There does 
exist, however, one paper which presents calculations of fin-line for both propagating 
and evanescent modes [5J. In order to adapt our analysis to fin-line, we need only 
change the dyadic Green’s function, Z i} . This was done, and a dispersion curve was 
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calculated, which is shown in Figure 6. Clearly, satisfactory agreement with the 
results in Reference 5 has been achieved. Since the microstrip calculation is quite 
similar to the fin-line calculation, it is believed that the microstrip mode calculations 
should also be reasonable. 

Once the modes have been found and verified, the next step is to calculate the 
characteristic impedance of the uniform line. This is done for two reasons. First, it 
gives added confirmation that the dominant mode calculation is accurate. Second, it 
gives a first-order appro xim ation to the input impedance at a discontinuity. This result 
is useful as a comparison for results obtained with the mode matching technique. Since 
the characteristic impedance is defined only for TEM lines, one must be careful to 
choose a definition of characteristic impedance which is useful experimentally. The 
definition most commonly chosen is the power-current definition 


Z • 

_ » 

(5a) 

II 


P = V 2 Re // E XH’ dx dy 

(5.b) 

S 

I = / J z (x)dx 

(Sx) 


An alternative definition is the voltage-current definition 

Z„ = j (6a) 

o 

V =— J E y (yUy z =0 (6.b) 

Calculations were performed with both of these definitions, and are plotted in Figure 7. 
They correspond very well with data already published in Reference 4. 




Figure 6. Dispersion curve of fin-line. Calculated with one basis function and 101 
spectral terms. 





(Sl/VIO) 30NVQ3cllM I 


Figure 7. Characteristic impedance of microstrip using two definitions. Calculated 
with one basis function and 101 spectral terms. 
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13. Analysis of Discontinuities with Mode Matching 

Among the techniques for studying waveguide discontinuities, mode matching 
seems one of the most promising [6J. This technique is suitable for calculating the 
characteristics of abrupt discontinuities, such as those in Figures 2 and 4. In addition, 
this method is easily adaptable to a cascade of abrupt discontinuities by using 
generalized S-parameters. Finally, it may be possible to model a smooth taper as a 
cascade of step discontinuities. Hence, the mode matching technique is capable of 
solving a wide variety of problems. 

The mode matching technique begins by using the modes of the uniform 
microstrip calculated in the previous section. This involves finding the transverse fields 
for each mode and in each of the two waveguides, 


«io (n,y) = e xia (n,y) x + e yia (n, y)y 

(7a) 

h^Cn.y) = h xia (n,y) x + h yia (n, y)y 

(7.b) 


where i denotes the mode number, and a denotes that we are referring to the input 
waveguide. Similar expressions hold for the output waveguide. At the point of the 
discontinuity, these transverse fields must be continuous. Enforcing this boundary 
condition, one obtains 

oo oo 

(l+p)a,e al + £<!,•«„■ = £ bje bj (8a) 

;=2 j=i 

___ oo __ oo 

(l p) 0. ih a 1 £ a i hai ~ 2 b j^bj (8.b) 

i=2 j—\ 

where p is the reflection coefficient of the dominant mode, and a, and bj are the 
amplitudes of the modes in waveguides a and b, respectively. At this point, one takes 
appropriate inner products of Equation (8) and forms a matrix equation for p and for 
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the mode coefficients a, and bj. Once these have been found, equivalent circuit 
parameters, Z„ and Y „ , are calculated as shown in Figure 8. The circuit has now been 
characterized completely. 

When using the above method, it is necessary to find an alternate method of 
calculating the equivalent circuit parameters. This gives approximate values for Z„ 
and Y n , against which one can compare the results derived from the mode matching 
technique. The only alternate method in the literature which treats this problem uses a 
static technique, which is valid only for lower frequencies. An example of this 
technique is shown in Reference 7. Therefore, a lower frequency case was chosen for 
the initial study, in order to be able to compare the results with other sources. If the 
validity of the mode matching technique can be demonstrated at lower frequencies, one 
may assume the technique is valid also at higher f requencies. 

By calculating a low frequency case, one is able to compare the results for Z„ 
and Y„ to calculations made with other methods. In this case, Z n is approximated by 
the normalized characteristic impedance of waveguide 2, relative to waveguide 1. The 
power-current definition of characteristic impedance is used, as shown in Equation (5). 
Furthermore, one can get an approximation to T„ by referring to the junction 
capacitances calculated in Figure 8. 

A sample case of a microstrip discontinuity calculation is shown in Figure 8- 
The equivalent circuit parameters, Z n and Y „ , are plotted as a function of the number 
of waveguide modes used. The convergence of the normalized input impedance, Z n , to 
its correct value has been demonstrated. The input admittance Y n , however, seems to be 
dependent upon the number of modes used in each guide. It is necessary to get Y n to 
converge in order to obtain the junction capacitance. It is likely that improved 
convergence of Y n can be obtained by improving the efficiency of the inner product 
calculation, thus allowing the use of more modes in the mode matching technique. 
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Figure 8. Variation of normalized input impedance and normalized susceptance with 
the number of waveguide modes. Calculated with one basis function and 
101 spectral terms. It is expected that Z n =0.985 and Y „ =128X10' . 
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There are many areas in this technique where additional study is needed. First, 
the sensitivity of the final answer to the number of spectral terms, the number of basis 
functions, and the number of modes used must be checked. Next, it is necessary to 
eliminate the dependency on the number of spectral terms used, by implementing an 
asymptotic form for the calculation of spectral terms of very large order. This change 
should reduce the computation time significantly. Finally, this technique should be 
applied to a variety of different structures such as those shown shown in Figures 1 and 
2, and other structures useful in transitions between printed circuits and rectangular 
waveguides. 

In summary a model for analyzing several types of transitions was described. 
The model involved the breakdown of transitions into a cascade of abrupt printed 

circuit discontinuities, whose S-parameters were determined with mode matching 

) 

techniques. The initial results presented in this report demonstrate the method to be a 
promising one, although further work is still needed. 
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2. CROSSTALK IN HIGH SPEED DIGITAL CIRCUITS AND INTERCONNECTIONS 
2.1. Introduction 

The present trend toward high density packaging and the implemention of 
high speed digital devices have led to increasing demands for the characterization of the 
digital circuits and interconnections from the point of view of cross-talk and distortion 
of signal pulses propagating in these circuits. Our objective in this effort has been to 
develop models for predicting the coupling in multiconductor transmission lines and 
estimating the signal distortion due to mismatch in these lines. The signal corruption 
engendered by the cross-talk and the reflection due to terminations can introduce false 
logic levels and undesired triggering of logic gates. Hence an accurate estimation of the 
signal distortion and the noise is important in designing high speed digital circuits. 

Since the introduction of electromagnetic couplers, several models for coupled 
lines have been proposed [8]. A frequency domain characterization [9] was performed 
using an integral equation technique to compute the capacitance and inductance matrix 
for coupled lines. Bryant and Weiss [10] established the relationship between the 
electrical and geometrical parameters of coupled microstrip pairs using the Green’s 
function approach. Krage and Haddad [ill, [ 12] determined the inductive and coupling 
coefficients as well as the directivity for various geometries of coupled lines. 
Hammerstad and Jensen [13] succeeded in implementing a model with errors less than 
that caused by physical tolerance and recently Kirshing and Jansen [14] reported 
reported frequency dependent expressions for microstrip pairs with excellent accuracy. 
Semiempirical equations for even and odd mode expressions for the design of such 
coupled structures were also derived [15]. Various authors extended the analysis to 
three-line couplers [16], [17] and, recently, Chan and Mittra [18] used a spectral- 
iterative technique based on the conjugate gradient method to compute the inductance 
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capacitance matrices of n -line microstrip structures. 

In the prediction of signal propagation, early contributors [l9]-{2l] used a weak 
coupling approximation to predict the transient response in coupled pairs . Paul [22], 
[23] derived matrix chain parameter identities for the analysis of multiconductor 
transmission lines in the frequency domain for sinusoidal excitation. Using a similar 
technique, Zysman and Johnson [24] derived an equivalent circuit for two coupled 
lines. Chang [25] implemented a circuit model to describe the transient response of an 
n -conductor system using a congruence transformer approach. Such a model was 
implemented on CAD programs [26] to analyze the cross-coupling response of 
interconnections in high speed digital circuits. The above and several other equivalent 
circuit techniques have been proposed to describe terminal variables of multiconductor 
systems. However, no closed-form expression relating the transients to the position 
along the lines has been made available. 

In the present study, we utilize the inductance and capacitance matrices of an 
n -line system terminated with passive linear impedances to derive the time- and 
position-dependent expressions for the transient response. The first step in our approach 
is to determine the modes of propagation of the network. Next, a Fourier integral 
analysis is applied to the system and the time-domain solutions are obtained. Finally, 
the accuracy of the method is evaluated by comparing the theoretical results with 
measurements. 

2.2. Formulation - Natural Modes 

The basic differential equations governing the propagation of voltages and 
currents in an n -line system can be written in a matrix form as 
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_ 0V =T 6I 

6* df 


(9) 


- iL =r^- 
d * » 


do) 


where V and I are column vectors of dimension n representing the line voltages and 
currents with [V], x = V, and, [I], j = 7, which are the voltage and current on the ith 
line respectively. L and C are square matrices of dimension n representing inductance 
and capacitance per unit length respectively, with [L], j if i = j and 

tU; =L>>ifi=y. is the self inductance per unit length of an isolated line and 
L, ( J n) is the mutual inductance per unit length between ith and jth conductors. 


Similarly, [C], , = C (s) + £ Cj m * (k =i ) if i—j, and [C], j - — C ,; m ) if i = j , where 

k =1 

C (l) is the self capacitance per unit length of a single isolated line and C-}™ \ the 
mutual capacitance between ith and jth lines. These self and mutual parameters can 
be measured using TDR techniques [27]. Since the lines are identical, L and C are 
symmetric matrices. Equations (9) and (10) can be combined to yield 


a"v = irfv 

Qx 2 Qt 2 


(11) 


jQfL =cl^L 

Qx 2 Qt 2 


( 12 ) 


In general, LC and CL are not equivalent; moreover, L and C being symmetric 
does not imply that LC or CL are symmetric. However it is known [23] that there 
exist square matrices E and H which premultiply both sides of (11) and (12) to yield 


6 2 EV 

d* 2 


= ELCE’ 1 


a 2 ^ 

# 2 


(13) 
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^2- = HCXH" 1 -^. (14) 

qx 2 a ? 2 

where ELCE -1 and HCLH -1 are diagonal matrices of order n . Since LC is the adjoint 
matrix of CL, it can be shown that 

ELCE" 1 = HCLH" 1 = A^ 2 (15) 

Therefore (13) and (14) can be rewritten as 


a’V. . 

-w “•sr 


(16) 


— 5~ — A m — 
Qx 2 Qr- 


(17) 


where 

V m * EV (18) 

l m = HI (19) 


Equations (16) and (17) give the relations between the modal voltage and 
current matrices, V m and 1^ is a diagonal eigenvalue matrix whose elements are 
given by {A tn ] i = l/v mi , where v mi is the velocity of propagation of the ith mode. The 
ith row of E and H are the voltage and current eigenvectors associated with the ith 
mode whose modal variables are given by [V m ], x and [I m ], j respectively. For the case 
where L and C are symmetric, there exist n distinct modes of propagation associated 
with n real eigenvalues and the set of linearly independent (real) eigenvectors forms a 
basis in n -dimensional space. In general, the entries of E and H depend on the entries 
of both L and C Otherwise, E and H are equivalent. The line variables can be 
recovered by using 
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I = E~ l V m 

(20) 

V = H-‘I m 

(21) 


If harmonic variation with time is assumed, the solutions to (16) and (17) have 
the form 




_js« 

e Vm A + 


+ J“* 

e v “b 


, jax 


( 22 ) 


l m (t jc ) = Z m 1 


jgx 

e v ' m a - 


+ j“* 
e Vm B 


gju> 


(23) 


± j“* jo* 

where € m is a diagonal matrix of order n with [6 “h = e ' m , w is the 

angular frequency and A and B are modal coefficient column vectors of order n 

associated with the forward and backward waves respectively. Making use of these 

equations in (9) and (10) for forward and backward waves separately, one gets 


Z m = A m - 1 ELH-‘ (24) 

Z m is the modal impedance matrix (dimension n ) which relates modal voltage 
and current waves. Analogously, a line or characteristic impedance matrix Z c can be 
defined by requiring V = Z C I for forward and backward waves separately. When 
combined with the above relations, one finds 

Z c = E-‘Z m H = E -1 A m -1 EL (25) 

Since in general, E and H are not equivalent, Z m and Z c contain off-diagonal 
elements which account for the interdependence of modal and state variables. 
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23. Transient Response 

The previous section presented the development of the expressions for the 
various matrices associated with a multiconductor transmission line system, the 
relationships between these matrices and the solution for a time-harmonic excitation. 
When the applied signals are arbitrary (nonsinusoidal, nonperiodic), the response will 
contain contributions from all frequencies. In this case, the coefficient vectors are 
functions of frequency and the time domain solutions for modal voltages and currents 
are obtained via integration in the frequency domain as follows: 


n (r A) = / 




jwx 


e b aW+ e "W 


e din 


(26) 




_JWX t Jfa»X 

e Vm A(oj)- e Vm bCw) 


e ■ / “ ? d h) 


(27) 


The above equations can be rearranged in terms of the forward and backward voltage 
waves V+ and V~ to read 


where 


V m (r^c) = V£t,Jc) + V^t^c) 


I m (t jc) = Z m 1 


V+0,*)- Vr(r,*) 




JUX 


v+(r^) = / e Vm A(o ))e^d0) 


JUX 


= f e Vro B(aj )e^d 


O) 


(28) 

(29) 


(30) 


(31) 
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+ jux 

The matrix products e m A(co) and e m B(a>) are the Fourier transforms of 
V + and V “ respectively, and, therefore, satisfy the relations 


JtdX 


e Vm A(oj) = ~~f V£t jc) e-J" dt 
2t t 


(32) 




e Vm BLu >) = yX ) e dt 

27T * 


(33) 


We then apply the boundary conditions at x =0 and x =1 (see Fig. 9). This 


yields 


Z s Z d -‘V s (t) = E-‘V m (r ,0) + ZsH-'IJt ,0) 


(34) 


0 = ET l V m (t J . ) - Z L H _1 I m (i JL ) 


(35) 


Z s and 2^ are the source and load impedance matrices respectively. They are 
constructed such that [Z S _1 ], 7 =(Z, U « ) ) -1 + £ (Z,) smJ )~‘ (& =0 if i=j, 

k =z 

[Z % - x \j=-(Z i ( j m) Y l if t=j and [Z^ 1 ],, =(Z, (i * 0" l + £ (Z^O' 1 Gk =i) if i=j, 

A =i 

— (Z, ( / m, ) _l if i=j. Z d is a diagonal matrix with [Z d ]„ =Z , isg ’. Z, (sg) and 
Z, (z?) are the resistances between ith line and ground at the source and load terminals 
respectively, whereas Z{j m ' and Z,( Zm ' are the resistances between ith and jth lines at 
the source and load terminals respectively. V s (t) is the column vector associated with 
the n arbitrary voltage sources at x = 0. A more general situation involving, in 
addition, voltage sources at the load terminals could be solved by using superposition 
and, hence, is not treated here. Taking the transforms of (34) and (35), and, using Eqs. 
(28M31), we obtain 
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Z a Z d _1 W s (a>) = E" 1 1 A(w) + B(oj) 


+ Z^E-'A* 


A(ct>) — B(«) 


(36) 


0 = E _1 


_J“L + J"L 

e Vm A(o>)+ e v ' m B(<u) 


— Z L L -1 E -1 A m 






e m A(oj) — e m B(w) 


(37) 


where W s is the Fourier transform of V s (t), defined by 


W s (o>) = V s (t) e~ }u * dt 

2ir J 


(38) 


Equations (36) and (37) can be rearranged and combined to give 


A(ai) = 


jli>l jti)l 

i n -r s e Vm r L e Vm 


-i 


TWiw) 


(39) 


-M M 

bc&») = - e Vm r L e Vm A(«j) 


(40) 


where l n designates the unit matrix of order n. The matrices r s , T L and T are 
respectively the source reflection coefficient, the load reflection coefficient and the 
transmission coefficient defined by 


r s = 1. + EZ s L -1 E -1 A m 


-i 


l n - EZ s L-‘E- 1 A m 


(41) 


r, = 


l n + EX L L~ 1 'E~ i A m 


-l 


l n - EZ L L -1 E _1 A m 


(42) 


T = l n + EZ 3 L -i E- 1 A m 


-l 


EZ s Z d 


-l 


(43) 


As can be anticipated, if the lines are terminated at the load end with a 
network equivalent to the characteristic impedance of the system 
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(Z]_ = Z c = E -1 A m _1 EL. See Eqs. (24) and (25)), then = 0 and no reflection occurs. 
This is the fa mili ar result of the single line case. Such a situation can also be achieved 
with the source impedance network to eliminate re-reflections. In the general case 
where reflections occur at both ends, T s and T L are not diagonal, and Eqs. (39) and (40) 
cannot be simplified further. However, one can rewrite A(oj) in a geometric seriest as 



M jml 

Kill jhll 

2 

A(ti>) = 

i n + r s e Vm r L e Vm + 

r s e Vm r L e Vm 

+ ••• 


TWiw) 


(44) 


A(<u) = 


k =0 


jtal fed 


in p n m 


r s e m r L e 


TW,(oj) 


Making use of this relation, V + and V m can be obtained using the inverse 
transform relations, (30) and (31) 


V+(r i x) = S( t-— )* 


i n + r s 8(t — — ) * r L s(t — — ) + 


= 8(t — — ) * S 

V 

m k =0 


r s 8(t — — )* r L s(t — — ) 


*(*) 


* TVit) 


TV,(t) (45) 


jc ) = - 8(t-— ) * r L 8(t — — ) * 
v m v m 


k =0 


r s 8(t — — ) * r L s(t — — ) 

% V m 


*(*) 


TVit) 


(46) 


In these expressions, 8(t ), 8(t ) and 8(t ) are diagonal matrices 

V m 

of order n with diagonal elements [8(t — — )L = 8(f — — ), [S(t — — )L = Sit — — ) 

V mi ^ mi 

and [8(t )]„ = 8 it ) respectively. Sit — r) being the unit impulse function 


t In a more rigorous manner.the necessary conditions for convergence of the geometric series must, a 
priori, be established. If the ratio for the series is diagonalizable, these conditions require that the magnitude of 
the associated eigenvalues be less than unity. In all cases however, convergence can be tested numerically be- 
fore applying the geometric series expansion. 
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occurring at t =t. The symbol * designates a convolution and the exponent ki*) 
indicates that the term within brackets is being convolved with itself k — 1 times, 
d/ (f )]*' ) = 1 and [ / it )] 1( ‘ ) = f it)). If we recall from transform theory 


Sit —u ) * 8(r — v ) = 8(r — [u +v ]) (47) 

8 it-u)* f it) = / it-u) (48) 


then the evaluations of (45) and (46) are reduced to a simple sequence of matrix 
multiplications in which the only time dependent factors are the delayed and 
attenuated versions of the applied signals represented by V s (t). For an arbitrary time 
t , the voltages and currents are given by 


Vit; c) = E~ 1 |v+(f^) + V^^) 


I(t^c) = L- 1 E- 1 A m 


V^^c)-V^(r^) 


(49) 

(50) 


In most practical cases, the source and load reflection cefficients are not very 
large; hence, only the first few terms of the infinite series need be retained in (45) and 
(46). However, in the event the mismatch is significant and the reflection coefficients 
are large, it becomes necessary to include more terms in the representation for V,J(r jc ) 
and V^it ). Equations (45) through (50) can be programmed on a computer to 
simulate the transient response of an n -line system for which the coupling parameters, 
the eigenvalues and the eigenvectors are known, and the terminations and the 
excitations are specified. 
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2 A. Application - Pulse Response Simulation 

Many digital applications involve situations in which one or several 
transmission lines are excited with a pulse. Waveform distortions and cross-coupled 
signals are then observed at every point along the lines. In order to implement a 
realistic model for such phenomena, a four-line stripline structure (see Fig. 10) was 
fabricated on glass-epoxy (dielectric constant = 4-5). The characteristic impedance of 
each line was designed to be 70 Cl. First the propagation characteristics were measured 
in order to determine the self and mutual parameters and set up inductance and 
capacitance matrices. The eigenvalues and eigenvectors associated with the matrix 
product LC were computed using an EMSL routine. Then, the lines were terminated as 
shown in Fig. 1 1 and the modal and state variables were computed using Eqs. (45X50). 
Finally the responses at x = 0 were plotted for each lin e and compared with the actual 
waveforms. As can be seen from Figs. 12 and 13, the. agreement is very good within 
experimental errors. For the shown configuration, 4 iterations were needed to attain 
sufficient agreement. Other simulations were performed for microstrip as well as 7-line 
structures. In all cases, satisfactory agreement was obtained between theory and 
experiment. 

The construction of the inductance and capacitance matrices requires the 
knowledge of the mutual parameters between any two conductors of the system. In 
practice, only the self and mutual parameters between adjacent lines need to be 
measured. Mutual coefficients between lines separated by another conductor can be 
approximated by extrapolation since they do not influence the response of the system 
significantly. From the geometry of these structures (microstrip and stripline) it can be 
anticipated that the mutual capacitance between two conductors decreases much faster 
than the mutual inductance as the number of conductors separating them and the 
spacing increase. This’ suggests that the cross-coupling in these multiconductor 
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Figure 12. Simulated waveforms of the voltage at the near end (x =0) for each line of 
the network of Fig. 11. Applied pulse magnitude: 2V, rise and fall time: 2 
ns, pulse width: 30 ns. 
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Figure 13. Experimental waveforms at x=0 observed at the near end of each line for 
the network of Fig. 11. 
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structures is predominantly inductive. 

The preceding simulations also used pulses with rise times in the order of 1 
nanosecond. The bandwidth associated with such speed is nearly 3 GHz. The variations 
of the line parameters were assumed to be negligible over this range. For much faster 
pulses however, the self and mutual coefficients become strong functions of frequency 
and the above analysis is no longer accurate. In this case the frequency dependence of 
the line parameters must be deter min ed before using a numerical approach. 

2.5. Conclusion 

Simple closed-forms expressions for the propagation of transients in parallel 
conductors have been derived. These expressions describe the voltage and current at 
any point of a multiconductor system for which the coupling parameters are known. 
The analysis assumes that the terminations are real and passive and that the line 
parameters are invariant with frequency or time. Various experiments using pulse 
excitation were performed on microstrip and stripline structures and the agreement 
was found to be excellent for rise times and fall times less than one nanosecond. 
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